Percolation transition in quantum Ising and rotor models with sub-Ohmic dissipation 
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We investigate the influence of sub-Ohmic dissipation on randomly diluted quantum Ising and 
rotor models. The dissipation causes the quantum dynamics of sufficiently large percolation clusters 
to freeze completely. As a result, the zero-temperature quantum phase transition across the lattice 
percolation threshold separates an unusual super-paramagnetic cluster phase from an inhomogeneous 
ferromagnetic phase. We determine the low-temperature thermodynamic behavior in both phases 
which is dominated by large frozen and slowly fluctuating percolation clusters. We relate our results 
to the smeared transition scenario for disordered quantum phase transitions, and we compare the 
cases of sub-Ohmic, Ohmic, and super-Ohmic dissipation. 
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I. INTRODUCTION 

The interplay between geometric, quantum, and ther- 
mal fluctuations in randomly diluted quantum many- 
particle systems leads to a host of unconventional low- 
tcmpcrature phenomena. These include the singular 
thermodynamic and transport properties in quantum 
Griffiths phases 1,2 as well as the exotic scaling behav- 
ior of the quantum phase transitions between different 
ground state phases. 3,4 Recent reviews of this topic can 
be found, e.g., in Refs. 5 and 6. 

An especially interesting situation arises if a quantum 
many-particle system is diluted beyond the percolation 
threshold p c of the underlying lattice (see, e.g., Ref. 7 
and references therein). Although the resulting percola- 
tion quantum phase transition is driven by the geometric 
fluctuations of the lattice, the quantum fluctuations lead 
to critical behavior different from that of classical per- 
colation. In the case of a diluted transverse-field Ising 
magnet, the transition displays exotic activated (expo- 
nential) dynamic scaling 8 similar to what is observed 
at infinite-randomness critical points. 3,4 The percolation 
transition of the quantum rotor model shows conven- 
tional scaling (at least in the particle-hole symmetric case 
where topological Berry phase terms are unimportant 9 ), 
but with critical exponents that differ from their classi- 
cal counterparts. 10,11 For site-diluted Hciscnberg quan- 
tum antiferromagnets, further modifications of the criti- 
cal behavior were attributed to uncompensated geomet- 
ric Berry phases. 12 ' 13 

In many realistic systems, the relevant degrees of 
freedom are coupled to an environment of "heat-bath" 
modes. The resulting dissipation can qualitatively 
change the low-energy properties of a quantum many- 
particle system. In particular, it has been shown that 
dissipation can further enhance the effects of randomness 
on quantum phase transitions. In generic random quan- 
tum Ising models, for instance, the presence of Ohmic 
dissipation completely destroys the sharp quantum phase 
transition by smearing 14-19 while it leads to infinite- 



randomness critical behavior in systems with continuous- 
symmetry order parameter. 20-22 Interestingly, supcr- 
Ohmic dissipation does not change the universality class 
of random quantum Ising models 17,19 but plays a major 
role in systems with continuous-symmetry order param- 
eter. 23 

It is therefore interesting to ask what arc the effects of 
dissipation on randomly diluted quantum many-particle 
systems close to the percolation threshold. It has re- 
cently been shown that Ohmic dissipation in a diluted 
quantum Ising model leads to an unusual percolation 
quantum phase transition 24 at which some observables 
show classical critical behavior while others are modified 
by quantum fluctuations. 

In the present paper, we focus on the influence of sub- 
Ohmic dissipation (which is qualitatively stronger than 
the more common Ohmic dissipation) on diluted quan- 
tum Ising models and quantum rotor models. When cou- 
pled to a sub-Ohmic bath, even a single quantum spin 
displays a nontrivial quantum phase transition from a 
fluctuating to a localized phase 25 whose properties have 
attracted considerable attention recently (see, e.g., Ref. 
26 and references therein). Accordingly, we find that 
the quantum dynamics of sufficiently large percolation 
clusters freezes completely as a result of the coupling to 
the sub-Ohmic bath, effectively turning them into classi- 
cal moments. The interplay between large frozen clusters 
and smaller dynamic clusters gives rise to unconventional 
properties of the percolation transition which we explore 
in detail. 

Our paper is organized as follows: In Sec. II, wc define 
our models and discuss their phase diagrams at a quali- 
tative level. Section III is devoted to a detailed analysis 
of the quantum rotor model in the large- N limit where 
all calculations can be performed explicitly. In Sec. IV, 
we go beyond the large-A limit and develop a general 
scaling approach. We conclude in Sec. V. 
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II. MODELS AND PHASE DIAGRAMS 



B. Classical percolation theory 



A. Diluted dissipative quantum Ising and rotor 
models 



We consider two models. The first model is a d- 
dimensional (d > 2) site-diluted transverse-field Ising 
model 8, 27-29 given by the Hamiltonian 



Hi = -J ^2 'I- 'h 17 ■ 17 j ~ h *^2 Vi<>i , 



(1) 



(i,3) 



a prototypical disordered quantum magnet. The Pauli 
matrices erf and af represent the spin components at 
site i, the exchange interaction J couples nearest neigh- 
bor sites, and the transverse field h x controls the quan- 
tum fluctuations. Dilution is introduced via the random 
variables rji which can take the values and 1 with prob- 
abilities p and 1 — p, respectively. We now couple each 
spin to a local heat bath of harmonic oscillators, 16,30 



H 



H I+ E 



,n4,n a i,n + ^ A i,n°fK> + 



,t 



(2) 

where ai. n (a|„) is the annihilation (creation) operator 
of the n-th oscillator coupled to spin i; v^ n is its natural 
frequency, and A, iW is the coupling constant. All baths 
have the same spectral function 



5{uj 



0=2™^-^^-"/"% (3) 



with a and uj c being the dimensionless dissipation 
strength and the cutoff energy, respectively. The expo- 
nent C characterizes the type of dissipation; we are mostly 
interested in the sub-Ohmic case < ( < 1. For compar- 
ison, wc will also consider the Ohmic (C = 1) and super- 
Ohmic cases (£ > 1). Experimentally, local dissipation 
(with various spectral densities) can be realized, e.g., in 
molecular magnets weakly coupled to nuclear spins 31,32 
or in magnetic nanoparticles in an insulating host. 33 

The second model is a site-diluted dissipative quantum 
rotor model which can be conveniently defined in terms 
of the effective Euclidean (imaginary time) action 10 

A = j dr ^ J TjiTlj 4>i (t) • (j>j (t) + ^ Vi-Adyn [4>i] 



<<:/} 



Aiy„[</>] = f ^I>c W K| C <£K) • 0(~Wn) • (4) 



Here, the random variables rji = 0, 1 again implement 
the site dilution, and oj n are bosonic Matsubara frequen- 
cies. The rotor at site i and imaginary time r is de- 
scribed by 4>i(r): a ^-component vector of length TV 1 / 2 . 
Its Fourier transform in imaginary time is denoted by 
<p((jj n ). The dynamic action .Adyn stems from integrating 
out the heat-bath modes, with the parameter a measur- 
ing the strength of the dissipation, and the exponent £ 
characterizing the type of the dissipation, as in the first 
model [sec Eq. (3)]. 



We now briefly summarize the results of percolation 
theory 34 to the extent necessary for our purposes. 

Consider a regular c?-dimcnsional lattice in which each 
site is removed at random with probability p. 35 For small 
p, the resulting diluted lattice is still connected in the 
sense that there is a cluster of connected nearest neighbor 
sites (called the percolating cluster) that spans the entire 
system. For large p, on the other hand, a percolating 
cluster does not exist. Instead, the lattice is made up of 
many isolated clusters consisting of just a few sites. 

In the thermodynamic limit of infinite system volume, 
the two regimes are separated by a sharp geometric phase 
transition at the percolation threshold p = p c . The be- 
havior of the lattice close to p c can be understood as 
a geometric critical phenomenon. The order parameter 
is the probability Poo of a site to belong to the infinite 
connected percolation cluster. It is obviously zero in the 
disconnected phase (p > p c ) and nonzero in the percolat- 
ing phase (p < Pc)- Close to p c , it varies as 



\P~Pc 



(P < Pc) 



(5) 



where j3 c is the order parameter critical exponent of clas- 
sical percolation. (We use a subscript c to distinguish 
quantities associated with the lattice percolation tran- 
sition from those of the quantum phase transitions dis- 
cussed below). In addition to the infinite cluster, we also 
need to characterize the finite clusters on both sides of 
the percolation threshold. Their typical size, the corre- 
lation or connectedness length £ c , diverges as 



£c ~ \P ~ Pc 



(G) 



with v c the correlation length exponent. The average 
mass S c (number of sites) of a finite cluster diverges with 
the susceptibility exponent j c according to 



S c ~ \p - Pc 



(7) 



The complete information about the percolation criti- 
cal behavior is contained in the cluster size distribution 
n s , i.e., the number of clusters with s sites excluding the 
infinite cluster (normalized by the total number of lat- 
tice sites). Close to the percolation threshold, it obeys 
the scaling form 



n s (p) = s T "} '[(p - p c )s ac 



(8) 



Here, r c and a c are critical exponents. The scaling func- 
tion f(x) is analytic for small x and has a single maxi- 
mum at some x max > 0. For large |ar|, it drops off rapidly 
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(x > 0), (9) 
(x<0), (10) 



where B\ and B^ are constants of order unity. The clas- 
sical percolation exponents arc determined by r c and a c : 
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FIG. 1. (Color online) Schematic ground state phase diagram of the diluted dissipative quantum Ising model Eq. (2) for fixed 
values of £ < 1, ui c , and J. The three panels show three cuts through the three-dimensional parameter space of dilution p, 
transverse field h x , and dissipation strength a. (a) a-p phase diagram at a fixed transverse field h x with h x > hoc{a — 0) 
such that the dissipationless system is in the paramagnetic phase. This phase diagram also applies to the rotor model Eq. (4). 
(b) hx-p phase diagram at a fixed dissipation strength a. (c) h x -a phase diagram at fixed dilution p < p c . CSPM refers to 
the cluster super-paramagnetic phase, transition (i) denotes the smeared generic (field or dissipation-driven) quantum phase 
transition, and (ii) and (iii) denote the percolation quantum phase transitions in the two regimes with or without dynamic 
clusters, respectively. 



the correlation lengths exponent v c = (r c — 1)/ (da c ), the 
order parameter exponent /3 C = (t c — 2)/<x c , and the sus- 
ceptibility exponent j c = (3 — t c )/ct c . 

Right at the percolation threshold, the cluster size dis- 
tribution does not contain a characteristic scale, n s ~ 
s~ Tc , yielding a fractal critical percolation cluster of frac- 
tal dimension Df = dj(r c — 1). 



C. Phase diagrams 

Let us now discuss in a qualitative fashion the phase 
diagrams of the models introduced in Subsec. II A, begin- 
ning with the diluted dissipative quantum Ising model 
Eq. (2). If we fix the bath parameters £ and lu c and 
measure all energies in terms of the exchange interac- 
tion J, wc still need to explore the phases in the three- 
dimensional parameter space of transverse field h x , dis- 
sipation strength a and dilution p. A sketch of the 
phase diagram is shown in Fig. 1. For sufficiently large 
transverse field and/or sufficiently weak dissipation, the 
ground state is paramagnetic for all values of the dilu- 
tion p. This is the conventional paramagnetic phase that 
can be found for h x > h 00 (a) or, correspondingly, for 
ol < cxoo(h x ). Here, h 00 (a) is the transverse field at 
which the undiluted bulk system undergoes the transi- 
tion at fixed a while oioo(h x ) is its critical dissipation 
strength at fixed h x . 

The behavior for h x < h ao (a) [or a > oioo(hx)] depends 
on the dilution p. It is clear that magnetic long-range 
order is impossible for p > p Cl because the lattice consists 
of finite-size clusters that are completely decoupled from 
each other. Each of these clusters acts as an independent 
magnetic moment. For h x < h^ia) and p > p c , the 
system is thus in a cluster super-paramagnetic phase. 



Let us consider a single cluster of s sites in more de- 
tail. For small transverse fields, its low-energy physics 
is equivalent to that of a sub-Ohmic spin-boson model, 
i.e., a single effective Ising spin (whose moment is pro- 
portional to s) in an effective transverse- field h x (s) ~ 
h x e~~ Bs with B ~ hi( J/h x ) and coupled to a sub-Ohmic 
bath with an effective dissipation strength a s = so;. 8,24 
With increasing dissipation strength and/or decreasing 
transverse field, this sub-Ohmic spin-boson model un- 
dergoes a quantum phase transition from a fluctuat- 
ing to a localized (frozen) ground state. 25 This implies 
that sufficiently large percolation clusters are in the lo- 
calized phase, i.e., they behave as classical moments. 
The cluster super-paramagnetic phase thus consists of 
two regimes. If the transverse field is not too small, 
h\{a) < h x < /loo (&) [or if the dissipation is not too 
strong, a\{h x ) > a > Oioo(h x )], static and dynamic clus- 
ters coexist. Here, h\(a) is the critical field of a single 
spin in a bath of dissipation strength a while ai(h x ) is 
its critical dissipation strength in a given field h x . In 
contrast, for h x < h\(a) [or a > ai(h x )], all clusters are 
frozen, and the system behaves purely classically. 

Finally, for dilutions p < p c , there is an infinite- 
spanning percolation cluster that can support magnetic 
long-range order. Naively, one might expect that the 
critical transverse-field (at fixed dissipation strength a) 
decreases with dilution p because the spins are missing 
neighbors. However, in our case of sub-Ohmic dissipa- 
tion, rare vacancy-free spatial regions can undergo the 
quantum phase transition independently from the bulk 
system. As a consequence, the field-driven transition 
[transition (i) in Fig. 1] is smeared, 15,18 and the ordered 
phase extends all the way to the clean critical field /loo (a) 
for all p < p c . Analogous arguments apply to the critical 
dissipation strength at fixed transverse field h x . 
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The infinite percolation cluster coexists with a spec- 
trum of isolated finite-size clusters whose behavior de- 
pends on the transverse field and dissipation strength. 
Analogous to the super-paramagnetic phase discussed 
above, the ordered phase thus consists of two regimes. 
For hi(a) < h x < h^a) [or ai(h x ) > a > aoo(h x )], 
static (frozen) and dynamic clusters coexist with the 
long-range-ordered infinite cluster. For h x < h\{a) [or 
a > ai(h x )], all clusters are frozen, and the system be- 
haves classically. 

The phase diagram of the diluted quantum rotor model 
with sub-Ohmic dissipation (4) can be discussed along 
the same lines. After fixing the bath parameters C and w c 
and measuring all energies in terms of the exchange inter- 
action J, we are left with two parameters, the dilution p 
and the dissipation strength a. The zero-temperature be- 
havior of a single quantum rotor coupled to a sub-Ohmic 
bath is analogous to that of the corresponding quantum 
Ising spin. With increasing dissipation strength, the ro- 
tor undergoes a quantum phase transition from a fluc- 
tuating to a localized ground state. This follows, for 
instance, from mapping 36 the sub-Ohmic quantum ro- 
tor model onto a one-dimensional classical Hciscnberg 
chain with an interaction that falls off more slowly than 
1/r 2 . This model is known to have an ordered phase 
for sufficiently strong interactions. 37 As a result, all the 
arguments used above to discuss the phase diagram of 
the diluted sub-Ohmic transverse-field Ising model carry 
over to the rotor model Eq. (4). The a— p phase diagram 
of the rotor model thus agrees with the phase diagram 
shown in Fig. 1(a). 

In the following sections, we investigate the percolation 
quantum phase transitions of the models Eqs. (2) and 
(4), i.e., the transitions occurring when the dilution p is 
tuned through the lattice percolation threshold p c . These 
transitions are marked in Fig. 1 by (ii) and (iii). 



III. DILUTED QUANTUM ROTOR MODEL IN 
THE LARGE- iV LIMIT 



In this section, we focus on the diluted dissipative 
quantum rotor model in the large- N limit of an infinite 
number of order-parameter components. In this limit, 
the problem turns into a self-consistent Gaussian model. 
Consequently, all calculations can be performed explic- 
itly. 



A. Single percolation cluster 

We begin by considering a single percolation cluster 
of s sites. For a > a^,, this cluster is locally in the 
ordered phase. Following Refs. 38 and 39, it can therefore 
be described as a single large-TV rotor with moment s 
coupled to a sub-Ohmic dissipative bath of strength a s = 



so.. Its effective action is given by 



a cS = tY^ 



1 r 



ip(u n )r n ip(-u) n ) - sH z (oj n )ip(-uj n ) 



( u ) 

where r„ = e + sacj* *|w n |^, ip represents one rotor com- 
ponent and H z is an external field conjugate to the order 
parameter. 

In the large- N limit, the renormalizcd distance e from 
criticality of the cluster is fixed by the large- N (spher- 
ical) constraint (\iP(t)\ 2 ) = 1. In terms of the Fourier 
transform, ip{uJ n ) defined by 



%p(r) =Ty^jj(u tl ) exp [-kj„t] , 



(12) 



the large- N constraint for a constant field H z becomes 

2 



i 



l — c i 

SOLUlc W n 



K 



e 



1. 



(13) 



Solving this equation gives the renormalized distance 
from criticality e as a function of the cluster size s. 

At zero temperature and field, the sum over the Mat- 
subara frequencies turns into an integration, and the con- 
straint equation reads 



duj 



1 



6q + sauic 



= 1. 



(14) 



(We denote the renormalizcd distance from criticality at 
zero temperature and field by eo.) The critical size s c 
above which the cluster freezes can be found by setting 
eo = and performing the integral (14). This gives 



= 1/ [na(l - 0} 



(15) 



As we are interested in the critical behavior of the clus- 
ters, we now solve the constraint equation for cluster 
sizes close to the critical one, s c — s <C s c . This can 
be accomplished by subtracting the constraints at s and 
s c from each other. We need to distinguish two cases: 
1/2 < £ < 1 and ( < 1/2. In the first case, the resulting 
integral can be easily evaluated after moving the cut-off 
uj c to infinity. This gives 

e = as c [-Csin(7r/C)a(s c - s)fl^u c (for C > 1/2). 

(16) 

In the second case, £ < 1/2, we can evaluate Eq. (14) via 
a straight Taylor expansion in (s c — s). This results in 

e =a 2 s c ^(l-2C)(s c -sK (for (< 1/2). (17) 

It will be useful to rewrite Eqs. (16) and (17) in a more 
compact manner: 



e ( S ) = [A;(l-s/ Sc )r /(1 - x W 



(18) 



where A c = -(as c ) 1 / c C sin(7r/C) for C > 1/2, and A ( 
(as c ) 2 7r(l - 2C) for C < 1/2, and x = max{l/2,C}. 
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In order to compute thermodynamic quantities, we will 
also need the value of e(s) at non zero temperature. The 
constraint equation for small but nonzero temperature 
can be obtained by keeping the cj„ = term in the fre- 
quency sum of Eq. (13) discrete, while representing all 
other modes in terms of an w-integral. This gives 



T 
e 



cIlo- 



1 



1. 



(19) 



Solving this equation for asymptotically low tempera- 
tures results in the following behaviors. For clusters 
larger than the critical size, s > s c , e vanishes linearly 
with T via e = Tsj (s— s c ). Clusters of exactly the critical 
size have e — A x uj\~ x T x . For smaller clusters (s < s c ), 
low temperatures only lead to a small correction of the 
zero-temperature behavior eo- Writing e(T) = eo + ST, 
we obtain S = [s/(s c — s)][x/(l — as)]. Clusters with sizes 
close to the critical one show a crossover from the off- 
critical to the critical regime with increasing T. For 
s < s r . this means 



e(T) 



e a (l + 5T/e ) (for e > e T ), 



(otherwise) 



(20) 



with e T = A x ul~ x T x . 

The constraint equation at zero temperature but in a 
nonzero ordering field H z can be solved analogously. 39 
For asymptotically small fields, we find c(H z ) = 
sH z [s/ (s — Sc)] 1 / 2 in the case of clusters of size s > s c . At 
the critical size, e{H z ) = [^^(s,^) 21 ] 1 /' 1 ^', and 
for s < s c we obtain c{H z ) = eo + S(sH z ) 2 /eo- Larger 
fields lead to a crossover from the off-critical to the crit- 
ical regime. For s < s c , it reads 



<H Z ) 



e [l + S {sH z /e ) 2 } (for e > e H J, 



CR-_ 



(otherwise), 



(21) 



with e Hz = [A^ujl- x (sH z ) 2x ]y ( . 1+x \ 

Observables of a single cluster can now be determined 
by taking the appropriate derivatives of the free energy 
F c i = -T]n(Z) with 



(22) 



where 



T 



i-C 



cxp 



T sH z (uj n )sH z {-u n ) 



i-Ci \c 



(23) 

The dynamical (Matsubara) susceptibility and magneti- 
zation are then given by 



l — C i if 



and 



m c i{u> n ) = T 



s 2 H z {uj n ) 
l— C i \c 



(24) 



(25) 



respectively, where e is given by the solution of constraint 
equation discussed above. (Note that the contribution of 
a cluster of size s to the uniform susceptibility is pro- 
portional to s 2 ). Therefore, in the above two limiting 
cases, we can write the uniform and static susceptibility 
of a cluster of size s < s c as a function of temperature as 
follows 



Xcl (T)^s 2 /e(T). 



(26) 



Large clusters (s > s c ) behave classically, \ci ~ s(s — 
s c )/T, at low-temperatures. Finally, for the critical ones 
Xd « s 2 /e T - 

In order to calculate the retarded susceptibility \ c i(ui), 
we need to analytically continue the Matsubara suscep- 
tibility by performing a Wick rotation to real frequency, 
iuj n — , lj + iO. The resulting dynamical susceptibility 
reads 

s 2 

e + auj c ^ [cos(7r£/2) — i sin(7rC/2)sgn(w)] 

(27) 

Using Eq. (21), the single cluster magnetization in a 
small ordering constant field H z is given by 

u^j H ^ 2 / e o (for eo » e H J, 
m ci = XdH z w < 2 (28) 
I H z s /£h z (otherwise). 

Thermal properties (at zero field) can be computed 
by using the "remarkable formulas" derived by Ford et 
al., which express the free energy (the internal energy) 
of a quantum oscillator in a heat bath in terms of its 
susceptibility and the free energy (internal energy) of the 
free oscillator. For our model, they read, respectively 



F c i = -fx 
and 

U c i = — M 



i r°° 

- / du]F f (uj,T)lm 
n Jo 



du)Uf(w,T)Ti 



— lnxci( w ) 
auj 



— In Xd{u) 



(29) 



(30) 



Here, F f (u,T) = Tln[2sinh(w/(2T))] and U f (oj,T) = 
(w/2) coth(oj/(2T)). The extra \i terms stem from the 
Lagrange multiplier enforcing the large- N constraint. 39 

The entropy S c i = {U c i — F c {)/T can be calculated 
simply by inserting Eq. (27) into Eqs. (29) and (30) and 
computing the resulting integral. For the dynamical clus- 
ters (s < s c ), the low-temperature entropy behaves as 



Bqolslu}. ^ 



(31) 



where is a ^-dependent constant. At higher temper- 
atures (greater than T* ~ e^/^uJc 1 ^), the entropy be- 
comes weakly dependent on T. 41 

In the low-T limit, the specific heat C c i = T(dS c i/dT) 
thus behaves as 



B c (asul-(—. 

Co 



(32) 
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B. Complete system 

After discussing the behavior of a single percolation- 
cluster, we now turn to the full diluted lattice model. 
The low-energy density of states of the dynamic clusters 
Pdy(e) = J2 s <s n sH £ ~ e oi s )) i g obtained combining the 
single-cluster result Eq. (18) with the cluster-size distri- 
bution Eq. (8), yielding 

Pdy (e) = A- 1 - 1) !Mflf£ ( ±) {1 ~ 2X)/X , (33 ) 



UJ, 



where s(e) is the size of a cluster with renormalizcd dis- 
tance e from criticality [which can be obtained inverting 
Eq. (18)]. Notice that pd y shows no dependence on e in 
the case £ < 1/2. In particular, it does not diverge with 
e — > 0, in contrast to the case £ > 1/2. 

We now discuss the physics at the percolation transi- 
tion, starting with the total magnetization m. We have to 
distinguish the contributions m ( i y from dynamical clus- 
ters, m s t from frozen finite-size clusters, and moo from 
the infinite percolation cluster, if any. For zero order- 
ing field H z , mdy vanishes, because the dynamic clusters 
fluctuate between up and down. The frozen finite-size 
clusters individually have a non-zero magnetization, but 
it sums up to zero (m s t = 0), because they do not align 
coherently for H z = 0. Hence, the only coherent con- 
tribution to the total magnetization is m^. Since the 
infinite cluster is long-range ordered for small transverse 
field h x < hooi^a), its magnetization is proportional to 
the number of sites in the infinite cluster, giving 



m = m c 



Poc(p) 



\p-Pc\ Pc (for p<p c ), 
(forp>]5 c ). 



(34) 



The magnetization critical exponent ft is therefore given 
by its classical lattice percolation value ft c . In response 
to an infinitesimally small ordering field H z , the frozen 
finite-size clusters align at zero temperature, leading to a 
jump in m(H z ) at H z = 0. The magnitude of the jump is 



given by m st = 
m st ~ (1 - p c 



E 



s>s c 



At the percolation threshold, 



and it vanishes exponentially for 
both p — > and p — > 1. The total magnetization in 
an infinitesimal field (given by m^ + m s t) is analytic at 
p = Pc and only clusters with sizes below s c are not 
polarized. 

To estimate the contribution m dy of the dynamic clus- 
ters, we integrate the magnetization of a single cluster 
Eq. (28) over the DOS given in Eq. (33). For C > 1/2, 
we find that 

3(1-0/(1+0 



mdy = Cqu s 



H- 



(35) 



where n s is the density of critical clusters, and C, 

V K/(1 '°C/(2C 



mdy 



At 



- 1). For £ < 1/2, the integration gives 

SrHz 



1 + ln 



A) 



(A ( u, c s*H*y/3 



(36) 




FIG. 2. (Color online) The magnetization as a function of di- 
lution p for different ordering field H z at absolute zero. The 
solid line is the magnetization at H z = (The contribution of 
the infinite cluster only). The dashed line is for an infinites- 
imal field and the remaining ones represents stronger fields. 
Insets dysplay the histerisis curves in the (i) ordered and (ii) 
disordered phases. 



where #o is a cut-off energy. 

Because the three contributions to the magnetization 
have different ficld-dcpcndcncc, the system shows un- 
conventional hysteresis effects. The infinite cluster has 
a regular hysteresis loop (for p < p c ), the finite-size 
frozen clusters do not show hysteresis, but they con- 
tribute jumps in m(H z ) at H z = 0, and the dynamic 
clusters contribute a continuous but singular term (see 
Fig. 2). 

The low-temperature susceptibility is dominated by 
the contribution \st of the static clusters, with each one 
adding a Curie term of the form s(s — s c )/T. Summing 
over all static clusters, close to the percolation threshold, 
we find that 



Xst 



s>s c 



s{s - S c ) 



lf\P 



Pc 



(37) 



For p —> and p 1, the prcfactor of the Curie term 
vanishes exponentially The infinite cluster contribution 
Xoo remains finite (per site) for T — > 0, because the infi- 
nite cluster is in the ordered phase. 

To determine the contribution Xdy of the dynami- 
cal clusters, we integrate the single-cluster susceptibil- 
ity Eq. (26) over the low-energy DOS in Eq. (33). For 
£ > 1/2, this gives 



Xdy 



l-2f 



with C' c = A7^[£/(2£ - 1)]. For £ < 1/2, we find 



Xdy 



A -l n *c S C 



1 + ln 



7(1 



(A cWc T)i/2 



(38) 



(39) 



The retarded susceptibility of the fluctuating clusters 
can be obtained by integrating the single-cluster suscep- 
tibility Eq. (27) over the distribution Eq. (33), this leads 
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to 



lmxdy{UJ) = D c 



sgn( W ), (40) 



with Dq = 1 



(I-l) 7 rsm(0(I-2))/[sin(f)( 7 r(l-C))i- 2 ]. 
We notice that Im Xd y has no w-dependence for £ < 1/2. 

Finally, we consider the heat capacity. The dynami- 
cal cluster contribution can be obtained by summing the 
single-cluster heat capacity Eq. (32) over pd y (e), yield- 
ing C dy ~ n Sc s c (T/w c ) 1_C for £ > 1/2 and C dy ~ 
n Sc s c {T/u c ) C for £ < 1/2. 



IV. BEYOND THE LARGE- AT LIMIT: SCALING 
APPROACH 

In the last subsection, we have studied the percolation 
quantum phase transition of the diluted sub-Ohmic rotor 
model Eq. (4) in the large- ./V limit. Let us now go beyond 
the laige-N limit and consider the rotor model with a 
finite number of components as well as the quantum Ising 
model Eq. (2). 

We begin by analyzing a single percolation cluster of 
s sites. For strong dissipation a > (or weak fluctu- 
ations h x < hoc), this cluster can be treated as a com- 
pact object that fluctuates in (imaginary) time only. As 
pointed out in Sec. II C, in the presence of sub-Ohmic 
dissipation, such a cluster undergoes a continuous quan- 
tum phase transition from a fluctuating to a localized 
phase as a function of increasing dissipation strength or, 
cquivalently, cluster size s. 

Even though the critical behavior of this quantum 
phase transition is not exactly solvable, we can still write 
down a scaling description of the cluster free energy 

F cl {r,H z ,T) = b- 1 F cl {rb 1 '^ z '\H z b- ! »,Tb) (41) 

where r = a s — a c = (s — s c )a is the distance from 
criticality, b is an arbitrary scale factor, and v s z s and y s 
are the critical exponents of the single-cluster quantum 
phase transition. (We use a subscript s to distinguish 
the single-cluster exponents from those associated with 
the percolation quantum phase transition of the diluted 
lattice.) 

Normally, one would expect the two exponents v s z s 
and y s to be independent. However, because the sub- 
Ohmic damping corresponds to a long-range interaction 
in time, the exponent 77 takes the mean-field value 2 — £ 
for all £. 42-44 This also fixes the exponent y s in Eq. (41) 
to be y s = (l + C)/2. Thus, there is only one independent 
exponent in addition to £; in the following we choose the 
susceptibility exponent y s . This implies, via the usual 
scaling relations, that the correlation time exponent is 
given by v s z s = 7 s /£. 

The values of the cluster exponents in the large- N case 
of Sec. Ill are given by 7 S = £/(l — £) and v s z s = 1/(1 — 
£). In the general case of finite- N rotors and for the 



quantum Ising model, they can be found numerically 
Notice the scaling form of the free energy Eq. (41) applies 
to bath exponents £ > 1/2. For £ < 1/2, the single- 
cluster critical behavior is mean-field-like. 

The behavior of single-cluster observables close to the 
(single-cluster) quantum critical point can now be ob- 
tained by taking the appropriate derivatives of the free 
energy Eq. (41). For example, the static magnetic sus- 
ceptibility at T = and H z = behaves as 



X(r,u = 0) ~ r~ 



(42) 



Using this result, we can derive a generalization of the 
probability distribution Pd y {t) of the inverse static sus- 
ceptibilities e = x -1 . We find 



Pdy{e) 



dsn s S [e — c(s c — s) 7s ] 



c (l-7.)/7. 



(43) 

right at the percolation threshold. In the large- N limit, 
7 S = £/(l - C) implying p dy (e) - eC 1 - 2 ?)/? in agreement 
with the explicit result in Eq. (33). 

Let us now discuss how the properties of the percola- 
tion quantum phase transition in the general case differ 
from those obtained in the large- N limit in Sec. Ill B. We 
focus on the case £ > 1/2. If the single-cluster critical 
behavior is of mean- field type (£ < 1/2), the functional 
forms of the results in Sec. Ill B are not modified at all. 
The total magnetization is the sum of the magnetiza- 
tion TOqo of the infinite percolation cluster, m st stem- 
ming from the large (s > s c ) frozen percolation clus- 
ters, and m e i y provided by the dynamic clusters having 
s < s c . Both TOqo and m s t are completely independent of 
the single-cluster critical behavior. The behavior of the 
spontaneous (zero-field) magnetization across the perco- 
lation transition in the general case is thus identical to 
that of the large- TV limit [see Eq. (34) and Fig. 2]. In 
contrast, the magnetization-magnetic field curve of the 
dynamic clusters does depend on the value of 7,,. Inte- 
grating the single cluster-magnetization of all dynamic 
clusters [in analogy to Eq. (28)] gives 



m dy 



ff [l-C+2f/ 7 ,]/(l+C)_ 



(44) 



In the large- N limit, this recovers the result Eq. (35), as 
expected. 

The low-temperature susceptibility can be discussed 
along the same lines. The contributions Xoo 

and Xst do 

not depend on the single-cluster critical behavior. Inte- 
grating the single-cluster susceptibility over all dynamic 
clusters using (43) yields (at p = p c ) 



Xd y 



y(i — rs)C/7s 



(45) 



If we use the large- N value of 7^, we reproduce Eq. (38). 

The scaling ansatz Eq. (41) for the single-cluster free 
energy thus allows us to discuss the complete thermody- 
namics across the percolation quantum phase transition. 
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Dynamic quantities can be analyzed in the same man- 
ner. For example, the scaling form of the single-cluster 
dynamic susceptibility reads 

Xd(r,H z ,T,uj) = b 2 ^- 1 Xcl {rb 1 ' ( - v ' z '\H z b y %Tb,ujb) 

(46) 

The contribution of the fluctuating clusters to the low- 
temperature dynamic susceptibility can be found by in- 
tegrating the single-cluster contribution over the distri- 
bution Eq. (43). This leads to 

Im Xti » ~ M< 1 -T"Wt" sgnH . (47) 

In the large- TV limit this corresponds to Imxdyiu) ~ 
|w| 1-2 ^ sgn(w) in agreement with Eq. (40) for £ > 1/2. 

In summary, even though the critical behavior is not 
exactly solvable for finite--/V rotors and quantum Ising 
models, we can express the properties of the percolation 
quantum phase transition in terms of a single indepen- 
dent exponent of the single-cluster problem (which can 
be found, e.g., numerically). 

V. CONCLUSIONS 

We have investigated the effects of local sub-Ohmic 
dissipation on the quantum phase transition across the 
lattice percolation threshold of diluted quantum Ising 
and rotor models. Experimentally, such local dissipation 
(with various spectral densities) can be realized, e.g., in 
molecular magnets weakly coupled to nuclear spins 31 ' 32 
or in magnetic nanoparticles in an insulating host. 33 
Further potential applications include diluted two-level 
atoms in optical lattices coupled to an electromagnetic 
field, random arrays of tunneling impurities in crystalline 
solids or, in the future, large sets of coupled qubits in 
noisy environments. 

As even a single spin or rotor undergoes a localiza- 
tion quantum phase transition for sufficiently strong sub- 
Ohmic damping, the quantum dynamics of large perco- 
lation clusters in the diluted lattice freezes completely. 
The coexistence of these frozen clusters which effectively 
behave as classical magnetic moments and smaller fluc- 
tuating clusters, if any, leads to unusual properties of 
the percolation quantum phase transition. In this final 
section, we put our results into broader perspective. 

Let us compare the three different quantum phase tran- 
sitions separating the paramagnetic and ferromagnetic 
phases [transitions (i), (ii), and (iii) in Fig. 1]. The 
generic transition (i) occurs as a function of transverse 
field or dissipation strength for p < p c . This transi- 
tion is smeared by the mechanism of Rcf. 15 because 
rare vacancy-free spatial regions can undergo the quan- 
tum phase transition independently from the bulk sys- 
tem. For p < p c , these rare regions are weakly coupled 
leading to magnetic long-range order instead of a quan- 
tum Griffiths phase. 18,19 

In contrast, the percolation transitions (ii) and (iii) 
are not smeared but sharp. The reason is that different 



percolation clusters are completely decoupled for p > p c . 
Thus, even if some of these clusters have undergone the 
(localization) quantum phase transition and display local 
order, their local magnetizations do not align, leading to 
an incoherent contribution to the global magnetization. 
Deviations from a pure percolation scenario change this 
conclusion. If the interaction has long-range tails (even 
very weak ones), different frozen clusters will be coupled, 
and their magnetizations align coherently This leads to 
a smearing of the dilution-driven transition analogous to 
that of the transition (i). However, if the long-range tail 
of the interaction is weak, the effects of the smearing 
become important at the lowest energies only. What is 
the difference between the percolation transitions (ii) and 
(iii) in Fig. 1? If all percolation clusters are frozen [tran- 
sitions (iii)] low-temperature observables behave purely 
classically. If large frozen and smaller dynamic clusters 
coexist [transitions (ii)] quantum fluctuations contribute 
to the observables at the percolation transition. 

We now compare the case of sub-Ohmic dissipation 
considered here to the cases of Ohmic and super-Ohmic 
dissipation as well as the dissipationless case. To do so, 
we need to distinguish the quantum Ising model and the 
rotor model. 

The percolation transitions of the dissipationless and 
super-Ohmic rotor models display conventional critical 
behavior, but with critical exponents that differ from 
the classical percolation exponents. 38 (This holds for the 
particle-hole symmetric case in which complex Berry 
phase terms are absent from the action. 9 ) In the Ohmic 
rotor model, the percolation transition displays activated 
scaling as at infinite-randomness critical points. 38 

For the diluted quantum Ising model, the percolation 
transition displays activated scaling already in the dissi- 
pationless 8 and super-Ohmic cases. 19 In the presence of 
Ohmic dissipation, sufficiently large percolation clusters 
can undergo the localization transition independently 
from the bulk. The resulting percolation transition 24 
is similar to the one discussed in the present paper, it 
shows unusual properties due to an interplay of frozen 
and dynamic percolation clusters. 

All these results suggest that quantum phase tran- 
sitions across the lattice percolation threshold can 
be classified analogously to generic disordered phase 
transitions, 5,10 (provided the order parameter action does 
not contain complex terms). If a single finite-size per- 
colation cluster is below the lower critical dimension of 
the problem, it can not undergo a phase transition in- 
dependent of the bulk system. The resulting percolation 
transition displays conventional critical behavior (this is 
the case for the dissipationless and super-Ohmic rotor 
models). If a single finite-size cluster can undergo the 
transition by itself (i.e., it is above the lower critical di- 
mension of the problem), the resulting percolation tran- 
sition is unconventional with some observables behaving 
classically while others are influenced by quantum fluc- 
tuations. This scenario applies to the sub-Ohmic models 
studied in this paper as well as the Ohmic quantum Ising 
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model. Finally, if a single percolation cluster is right 
at the lower critical dimension (but does not undergo a 
phase transition), the percolation quantum phase transi- 
tion shows activated critical behavior. This scenario ap- 
plies to the dissipationless quantum Ising model as well 
as the Ohmic quantum rotor model. 
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